Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody (preprint)

The COVID-19 pandemic has led to over 760 million cases and 6.9 million deaths worldwide. To mitigate the loss of lives, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with susceptible variants. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response anti-viral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

How reliable are estimates of key parameters in viral dynamic models? (preprint)

Mathematical models of viral infection have been developed and fit to data to gain insight into disease pathogenesis for a number of agents including HIV, hepatitis C and B virus. However, for acute infections such as influenza and SARS-CoV-2, as well as for infections such as hepatitis C and B that can be acute or progress to being chronic, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the exponential phase of viral growth, i.e., when most transmission events occur. Missing data may make estimation of the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. Here, we evaluated the reliability of estimates of key model parameters when viral load data prior to the viral load peak is missing. We estimated the time from infection to peak viral load by fitting non-linear mixed models to a dataset with frequent viral RNA measurements, including pre-peak. We quantified the reliability of estimated infection times, key model parameters, and the time to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. We find a lack of data in the exponential growth phase underestimates the time to peak viral load by several days leading to a shorter predicted exponential growth phase. On the other hand, having an idea of the time of infection and fixing it, results in relatively good estimates of dynamical parameters even in the absence of early data.

An explanation for SARS-CoV-2 rebound after Paxlovid treatment (preprint)

In a fraction of SARS-CoV-2 infected individuals treated with the oral antiviral Paxlovid, the virus rebounds following treatment. The mechanism driving rebound is not understood. Here, we show that viral dynamic models based on the hypothesis that Paxlovid treatment near the time of symptom onset halts the depletion of target cells, but may not fully eliminate the virus, which can lead to viral rebound. We also show that the occurrence of viral rebound is sensitive to model parameters, and the time treatment is initiated, which may explain why only a fraction of individuals develop viral rebound. Finally, the models are used to test the therapeutic effects of two alternative treatment schemes. These findings also provide a possible explanation for rebounds following other antiviral treatments for SARS-CoV-2.

Risk of BA.5 infection in individuals exposed to prior SARS-CoV-2 variants (preprint)

The SARS-CoV-2 omicron BA.5 subvariant is progressively displacing earlier subvariants, BA.1 and BA.2, in many countries. One possible explanation is the ability of BA.5 to evade immune responses elicited by prior BA.1 and BA.2 infections. The impact of BA.1 infection on the risk of reinfection with BA.5 is a critical issue because adapted vaccines under current clinical development are based on BA.1. We used the national Portuguese COVID-19 registry to analyze the risk of BA.5 infection in individuals without a documented infection or previously infected during periods of distinct variants' predominance (Wuhan-Hu-1, alpha, delta, BA.1/BA.2). National predominance periods were established according to the national SARS-CoV-2 genetic surveillance data (when one variant represented >90% of the sample isolates). We found that prior SARS-CoV-2 infection reduced the risk for BA.5 infection. The protection effectiveness, related to the uninfected group, for a first infection with Wuhan-Hu-1 was 52.9% (95% CI, 51.9 - 53.9%), for Alpha 54.9% (51.2 - 58.3%), for Delta 62.3% (61.4 - 63.3%), and for BA.1/BA.2 80.0% (79.7 - 80.2%). The results ought to be interpreted in the context of breakthrough infections within a population with a very high vaccine coverage (>98% of the study population completed the primary vaccination series). In conclusion, infection with BA.1/BA.2 reduces the risk for breakthrough infections with BA.5 in a highly vaccinated population. This finding is critical to appraise the current epidemiological situation and the development of adapted vaccines.

A simple model of COVID-19 explains disease severity and the effect of treatments (preprint)

Considerable effort was made to better understand why some people suffer from severe COVID-19 while others remain asymptomatic. This has led to important clinical findings; people with severe COVID-19 generally experience persistently high levels of inflammation, slower viral load decay, display a dysregulated type-I interferon response, have less active natural killer cells and increased levels of neutrophil extracellular traps. How these findings are connected to the pathogenesis of COVID-19 remains unclear. We propose a mathematical model that sheds light on this issue. The model focuses on cells that trigger inflammation through molecular patterns: infected cells carrying pathogen-associated molecular patterns (PAMPs) and damaged cells producing damage-associated molecular patterns (DAMPs). The former signals the presence of pathogens while the latter signals danger such as hypoxia or the lack of nutrients. Analyses show that SARS-CoV-2 infections can lead to a self-perpetuating feedback loop between DAMP expressing cells and inflammation. It identifies the inability to quickly clear PAMPs and DAMPs as the main contributor to hyperinflammation. The model explains clinical findings and the conditional impact of treatments on disease severity. The simplicity of the model and its high level of consistency with clinical findings motivate its use for the formulation of new treatment strategies.

Emergence of SARS-CoV-2 Resistance with Monoclonal Antibody Therapy (preprint)

Resistance mutations to monoclonal antibody (mAb) therapy has been reported, but in the non-immunosuppressed population, it is unclear if in vivo emergence of SARS-CoV-2 resistance mutations alters either viral replication dynamics or therapeutic efficacy. In ACTIV-2/A5401, non-hospitalized participants with symptomatic SARS-CoV-2 infection were randomized to bamlanivimab (700mg or 7000mg) or placebo. Treatment-emergent resistance mutations were significantly more likely detected after bamlanivimab 700mg treatment than placebo (7% of 111 vs 0% of 112 participants, P=0.003). There were no treatment-emergent resistance mutations among the 48 participants who received bamlanivimab 7000mg. Participants with emerging mAb resistant virus had significantly higher pre-treatment nasopharyngeal and anterior nasal viral load. Intensive respiratory tract viral sampling revealed the dynamic nature of SARS-CoV-2 evolution, with evidence of rapid and sustained viral rebound after emergence of resistance mutations, and worsened symptom severity. Participants with emerging bamlanivimab resistance often accumulated additional polymorphisms found in current variants of concern/interest and associated with immune escape. These results highlight the potential for rapid emergence of resistance during mAb monotherapy treatment, resulting in prolonged high level respiratory tract viral loads and clinical worsening. Careful virologic assessment should be prioritized during the development and clinical implementation of antiviral treatments for COVID-19.

Kinetics of SARS-CoV-2 infection in the human upper and lower respiratory tracts and their relationship with infectiousness (preprint)

SARS-CoV-2 is a human pathogen that causes infection in both the upper respiratory tract (URT) and the lower respiratory tract (LRT). The viral kinetics of SARS-CoV-2 infection and how they relate to infectiousness and disease progression are not well understood. Here, we develop data-driven viral dynamic models of SARS-CoV-2 infection in both the URT and LRT. We fit the models to viral load data from patients with likely infection dates known, we estimated that infected individuals with a longer incubation period had lower rates of viral growth, took longer to reach peak viremia in the URT, and had higher chances of presymptomatic transmission. We then developed a model linking viral load to infectiousness. We found that to explain the substantial fraction of transmissions occurring presymptomatically, the infectiousness of a person should depend on a saturating function of the viral load, making the logarithm of the URT viral load a better surrogate of infectiousness than the viral load itself. Comparing the roles of target-cell limitation, the innate immune response, proliferation of target cells and spatial infection in the LRT, we found that spatial dissemination in the lungs is likely to be an important process in sustaining the prolonged high viral loads. Overall, our models provide a quantitative framework for predicting how SARS-CoV-2 within-host dynamics determine infectiousness and represent a step towards quantifying how viral load dynamics and the immune responses determine disease severity.